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^ New double ratios of exponential moments of different two-point functions, which are less sensitive 
Q"^ ■ to the heavy quark mass and to the continuum effects than the commonly used ratio of moments, 
are presented for a more refined analysis of the mass-splittings between the different heavy quarko- 
Ph! nia states. We show that at the c and b quark mass scales the 1/m-expansion does not converge 
Q^l for these quarkonia channels, while a connection of our mass and width formulae, with the poten- 
^ ' tial model ones is done. Using the present value of the QCD coupling as, we deduce the value: 
{asG^) = (7.5 ± 2.5) x 10"^ GeV^ of the gluon condensate from - M^^ and M^, - Mx, which 
we compare with the ones from different fits of the heavy and light quark channels. We also find 



that M^^(pi) — M^^(p3) is gouverned by the radiative corrections and gives as(1.3 GeV)=0.64 



+0.36 
-0.18 



I for 4 flavours, implying as{Mz) = 0.127 ± 0.011. Our predictions for the splittings of different 
heavy quarkonia states are summarized in Table 2, where, in particular, we find M-^ — M^^ ^ 
MeV implying the possible observation of the rjf, in the T-radiative decay. 
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The double ratio of moments 



QCD spectral sum rule (QSSR) a la SVZ (for a recent review, see e.g. 0]) has shown since 
15 years, its impressive ability for describing the complex phenomena of hadronic physics with 
the few universal "fundamental" parameters of the QCD Lagrangian (QCD coupling a^, quark 
masses and vacuum condensates built from the quarks and/or gluon fields), without waiting for a 
complete understanding of the confinement problem. In the example of the two-point correlator: 

nQ(g^) ^tjd'x e^^^ {0\TJq{x) {Jq{o))^ |0), (1) 

associated to the generic hadronic current: Jq{x) = QTQ{x) of the heavy Q-quark (F is a Dirac 
matrix which specifies the hadron quantum numbers), the SVZ-expansion reads: 



£(%Vi^/^)(0(/.)) 

D=0,2,... dimO=D ^Mq — 



no(g') -EE — , „ .D/2 — ' (2) 



where is an arbitrary scale that separates the long- and short- distance dynamics; C^"^^ are the 
Wilson coefficients calculable in perturbative QCD by means of Feynman diagrams techniques; 
(O) are the non-perturbative condensates of dimension D built from the quarks or/and gluon 
fields {D = corresponds to the case of the naive perturbative contribution). Owing to gauge in- 
variance, the lowest dimension condensates that can be formed are the D = 4 light quark mq{ipip) 
and gluon {agG"^) ones, where the former is fixed by the pion PCAC relation, whilst the latter is 
known to be (0.07 ± 0.01) GeV^ from more recent analysis of the light quark systems p. The 
validity of the SVZ-expansion has been understood formally, using renormalon techniques (absorp- 
tion of the IR renormalon ambiguity into the definitions of the condensates and absence of some 
extra 1/g^-terms not included in the OPE) ^] and/or by building renormalization-invariant 
combinations of the condensates (Appendix of and references therein). The SVZ expansion is 
phenomenologically confirmed from the unexpected accurate determination of the QCD coupling 
ttg and from a measurement of the condensates from semi-inclusive r-decays 0-0. 
The previous QCD information is transmitted to the data through the spectral function ImnQ(t) 
via the Kallen-Lehmann dispersion relation {global duality) obeyed by the hadronic correla- 
tors, which can be improved from the uses of different versions of the sum rules [|l|, ||9[|-[p 



In this paper, we shall use the simple duality ansatz parametrization: "one narrow resonance" + 
'^QCD continuum" , from a threshold tc, which gives a good description of the spectral integral in 
the sum rule analysis, as has been tested successfully in the light-quark channel from the e"'"e~ 
1 = 1 hadron data and in the heavy-quark ones from the e~^e~ ^ ip ot T data. We shall work 
with the relativistic version of the Laplace or exponential sum rules where the QCD expression 
konown to order as is given in terms of the pole mass m(p^ = rri^) [|ll|,|jl2[-|l^ 0: 



POO 

C{a,m'^) = / dt exp{-ta) - lmUQ{t) = Ani^AHiuj) 

Jim? TV 



1 + asaniuj) + — —bH{uj) 

36 



TZnia) = — ^ log Ch {a, m^) =Am^FH{uj) 
da 



(3) 



where: 



= 1/x = 4m^(T (4) 

is a dimensionless variable, while a = r (notation used in the literature) is the exponential Laplace 
sum rule variable; Fh, Ph and Qh are complete QCD Whittaker functions compiled in Ref. 



12|-|14|; H specifies the hadronic channel studied. In principle, the pair {(T,tc) is free external 



^For consistency, we shall work with the too- loop order ag expression of the pole mass [151. 



parameters in the analysis, so that the optimal result should be insensitive to their variations. 
Stability criteria, which are equivalent to the variational method, state that the best results should 
be obtained at the minimas or at the inflexion points in n or a, while stability in tc is useful to 
control the sensitivity of the result in the changes of tc-values. These stability criteria are satisfied 
in the heavy quark channels studied here, as the continuum effect is negligible and does not exceed 
1% of the ground state contribution 0, |12|, such that at the minimum in cr, one expects to a good 
approximation: 

mm^n{a) ~ M^. (5) 

Moreover, one can a posteriori check that, at the stability point, where we have an equilibrium 
between the continuum and the non-perturbative contributions, which are both small, the OPE 
is still convergent such that the SVZ-expansion makes sense. The previous approximation can be 
improved by working with the double ratio of moments : 



7^ 



7^ 



HH'{X 



H 



n 



H' 



Ml 



^0 



HH' 



y 



2 \HH' 



(6) 



provided that each ratio of moments stabilizes at about the same value of a, as in this case, there 
is a cancellation of the different leading terms such as the heavy quark mass (and its ambiguous 
definition used in some previous literatures), the negligible continuum effect (which is already 
small in the ratio of moments), and each leading QCD corrections. We shall limit ourselves here 
to the a^-correction for the perturbative contribution and to the leading order one in as for the 
gluon condensate effects. To the order we are working, the gluon condensate is well-defined as the 
ambiguity only comes from higher order terms in a^, which have, however, a smaller numerical 
effect than the one from the error of the phenomenological estimate of the condensate. 



Test of the 1/m-expansion 



For this purpose, we use the complete horrible results expressed in terms of the pole mass to 
order given by []T2[ and checked by various authors 0, which we expand with the help of the 
Mathematica program. We obtain for different channels the expressions given in Table 1. By 
comparing the complete and truncated series in 1/m, one can notice that, at the c and b mass 
scales, the convergence of the 1/m-expansion is quite bad due to the increases of the numerical 
coefficients with the power of 1/m and to the alternate signs of the 1/m series. 



Balmer-mass formula from the ratio of moments 



The Balmer formula derived from a non-relativistic approach (m — * oo) of the Schrodinger levels 
reads [|T^] (see also [|TH]-[^). for the Sf vector meson: 



'53 



2m 
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+ 0.23- 



TT 



masY 



(7) 



It is instructive to compare this result with the mass formula obtained from the ratio of moments 
within the 1/m-expansion. Using the different QCD corrections in Table 1, one obtains the mass 
formula at the minimum in cr of TZ: 



M53 



2m(l + 



a 



1 



6m 



cr (OsG ) 



12 



(8) 



*This method has also been used in |Iq] for studying the mass splittings of the heavy-light quark systems. 



Table 1: Expanded expressions of different QCD corrections in tfie case of the pole mass 
m^) known to order a^- 









Vector Sf 


TxAv 

av 
by 
Fy 
Pv 
Qv 


3 3/2 / 1 3 _|_ 45 ™2 525 ™3 , A 

4 ^ 0.040 + 1.952v^ - 1.539a; - ..) . 

1 I 3 _|_ 27 21 _|_ 
2x^ ' 2x'^ ' 8x 8 ' ■■■ 

1 + fx - fx^ + - ... 
-§0^ + 2.704x^/2 _ io.093a;5/2 + 52.93x^/2 - ... 

3 1 , 13 41 , 
2x2 4x ^ 8 4 ^ ■■■ 




S-waves splitting 


Ar 
Ar 
Ar 


1-^ + 1x3-... 

^^3/2 + 1.539x2 - 3.0258x-'^/2 _ ^^^iq^s ^ 26.307x^/2 ^ 




P-waves splittings 


A 01 

Aq 

^01 

a 

AOl 

A^ 

Ar 
Ar 


i 

-3.18x2(1 - 10.17X + 102.1x2 + ...) 

2 _|_ 5 55 ^ _|_ 
X 2 ^ ^ "T ... 

1.06x2(1 -9.5x + 81.1x2 - ...) 

1 + x2 - ^X^ + ... 

-0.1576x^/2 - 2.545x2 + 3.95x^/2 - ... 

6 I 31 89™ 1715^2 , 

4 8-^ gX-r... 




P- versus S-waves splittings 


AVS 
^0 

^vs 

a 

AVA 

^0 

^VA 

a 

Ar 
Ar 

^VT 

a 

AVT 


1 - X + 5x2 - 3Q^3 

-f^TTX - 0.336x^/2 + 4.244x2 + 7.458x^/2 - 42.017x3 - ... 

3 2 41 1 389 
a;2 X 4 4 

Ar 

-fv^TTX - 0.336x^/2 + 1.06x2 + 7.458x^/2 - 9.655x3... 

3 4 31 _|_ 167 
x2 X 4 2 
1 - X + 6x2 - ^^3 _^ 

-fr^ - 0.493x^/2 - 1.484x2 + 11.409x^/2 + 18.248x3 - ... 

3 10 1 579 16719 2 , 
x'^ X ^ 8 X -r ... 



In the case of the 6-quark, where we expect the static approximation more rehable,the minimum 
of this quantity is obtained to leading order at: 



V '^coul 



4:1710, 



0.85 GeV 



TT 



(9) 



where we have used for 5 flavours Q as{<j) — 0.32 ± 0.06 after evolving the value as{Mz) = 
0.118 ± 0.006 from LEP |2l|] and r-decay data The inclusion of the gluon condensate 



correction shifts the value of CTmm to: 



(10) 



These previous values of a confirm the more involved numerical analysis in and indicates the 
relevance of the gluon condensate in the analysis of the spectrum. By introducing the previous 
leading order expression of acoui into the sum rule, one obtains: 



2m,, 1 



TT 



27 
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27 
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3^2 
7r 



TT 
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where one can deduce by identification in the static limit {mi, oo) that the Coulombic effect is 
exactly the same in the two approaches. The apparent factor 7r/3 is due to the fact that we use 
here the approximate Schwinger interpolating formula for the two-point correlator . The gluon 
condensate coefficient is also about the same in the two approaches. This agreement indicates the 
consistency of the potential model and sum rule approach in the static limit, though a new extra 
al correction due to the v'^ (finite mass) terms in the free part appears here (for some derivations 
of the relativistic correction in the potential approach see ||2^, [^), and tends to reduce the 



coulombic interactions. On the other hand, at the 6-quark mass scale, the dominance of the gluon 
condensate contribution indicates that the 6-quark is not enough heavy for this system to be 
coulombic rendering the non-relativistic potential approach to be a crude approximation at this 
scale. 



Sf — Sq hyperfine and P — ^-wave splittings 



In the non-relativistic approach a la pO[, the hyperfine and S — P wave splittings are given to 
leading order by: 



M{Sf)-M{S^) ~ 2m 
M(Pf ) - M{Sl) - 2m 



32 



1 + 3.255 



TT 



2m%6 

{asG^} 



(12) 



where Cp = 4/3. Using the double ratio of moments and the QCD corrections given in Table 1, 
one obtains at acoui- 

M{S!)-M{Sl) ^ 



M{Pf) - M{Sf) 




M{Sf) 



47r 2 27r 2 27 , 



(13) 



^In the approximation we are working, the effect of the number of flavours enters only through the /3-function 
and therefore is not significant. 



where the corrections are, respectively, relativistic, Coulombic and non-perturbative. By compar- 
ing the sum rules in Eqs. (8) and (13), one can realize that the leading x or 1/cr-terms cancel in 
the hyperfine splitting, while the x-expansion is slowly convergent for the a^-term at the 6-mass. 
Comparing now this result with the one from the non-relativistic approach, it is interesting to 
notice that both approaches lead to the same Os-behaviour of the Coulombic and gluon conden- 
sate contributions. A one to one correspondence of each QCD corrections is not very conclusive, 
and needs an evaluation of the correlator at the next-next-to-leading order for a better control of 
the a^x terms. However, as the Coulombic potential is a fundamental aspect of QCD, we shall, 
however, expect that, after the resummation of the higher order terms in a^, the coefficient of 
the a^-term in the hyperfine splitting will be the same in the two alternative approaches. In the 
case of the S — P wave splitting, the sum of the a1 corrections agrees from the two methods, 
though one can also notice that the relativistic correction is larger than the Coulombic one. The 
discrepancy for the coefficients of the gluon condensate in the two approaches is more subtle and 
may reflect the difficulty of Bell-Bertlmann [|I^ to find a bridge between the field theory a la SVZ 
(fiavour-dependent conflning potential) and the potential models (flavour- independence). Resolv- 
ing the diffferent puzzles encountered during this comparison is outside the scope of the present 
paper. 



Leptonic width and quarkonia wave function 

Using the sum rule Ch and saturating it by the vector Sf state, we obtain, to a good approximation, 
the sum rule: 



(«eQ)^ 



^25mMvcT ^-3/2 



1 + --v/Trama, (a^G'^)ma^^'^ 

gV 9 ^ " / 



(14) 



where eg is the quark charge in units of e; 6m = My — 2m is the meson-quark mass gap. In the 
case of the 6-quark, we use [|T5[ 6m ~ 0.26 GeV, and the value of amin given in Eq. (10). Then: 



rx(s3)^e+e- ^ 1-2 keV , (15) 

in agreement with the data 1.3 keV. However, one should remark from Eq. (14), that the ag 
correction is huge and needs an evaluation of the higher order terms (the gluon condensate effect 
is negligible), while the exponential factor effect is large, such that one can reciprocally use the 
data on the width to fix either or/and the quark mass. Larger value of the heavy quark mass 
at the two- loop level (see e.g. [Q) corresponding to a negative value of 6m, would imply a smaller 
value of the leptonic width in disagreement with the data. 

In the non-relativistic approach, one can express the quarkonia leptonic width in terms of its wave 
function ^(0)q: 

We- = ^-^elimil (l - 4C^f ) , (16) 

where (see e.g. pO| ): 

fl-4C^"' 

TT 



167r|^(0)|^ (1-4Cf^) ~ 2(mC^a,)^ ^ 15 GeV^ (17) 
In our approah, one can deduce: 

167r|*(0)||(l-4C^^ 



1 ^2SmMvtT^-3/2 ^V 

72y/TT m 
18.3 GeV^ . 



8 , AtT 2\ 

1 + -y/'Kamas — ^i'^sC^ )ma ' 



Using the expression of <7coui, one can find that, to leading order, the two approaches give a similar 
behaviour for '^{0)q in as and in m and about the same value of this quantity, though, one should 
notice that in the present approach, the QCD coupling as is evaluated at the scale a as dictated by 
the renormalization group equation obeyed by the Laplace sum rule but not at the resonance 
mass! 



Gluon condensate from M 



The value of a, at which, the S-wave charmonium ratio of sum rules stabilize is: [p!2|: 



a - (0.9 ±0.1) GeV~ 



(19) 



Using the range of the charm quark pole mass to order as accuracy 0: ~ 1.2 — 1.5 GeV 
one can deduce the conservative value of x: 



1/x ~ 6.6 ± 1.8 



(20) 



The ratio of the mass squared of the vector V{Sf) and the pseudoscalar P{Sq) is controlled by 
the double ratio of moments given generically in Eq. (6), where the exact expressions of the 
corrections read: 



A 



VP 



0.995 



+0.001 
0.004 



A 



VP 



0.0233 



-0.009 
+0.011 



A 



VP 
G 



90 77+8-86 
ZJ. ( ( _io.23' 



(21) 



where each terms lead to be about 0.5, 2 and 7 % of the leading order one. One can understand 
from the approximate expressions in Table 1 that the leading x-corrections appearing in the ratio of 
moments cancel in the double ratio, while the remaining corrections are relatively small. However, 
the x-expansion is not convergent for the a^-term at the charm mass, which invalidates the use of 
the 1/m-expansion done in p8| in this channel. Using for 4 flavours [jl5|: as^a) ~ 0.48to;io, and 
the experimental data 



VP 



1.082, one can deduce the value of the gluon condensate: 

{asG^) ^ (0.10 ± 0.04) GeV^ (22) 

We have estimated the error due to higher order effects by replacing the coefficient of as with the 
one obtained from the effective Coulombic potential, which tends to reduce the estimate to 0.07 
GeV^. We have tested the convergence of the QCD series in a, by using the numerical estimate 
of the dimension-six gluon condensate g{fabcG"'G^G'^) contributions given in [|I4|. This effect is 



about 0.1% of the zeroth order term and does not influence the previous estimate in Eq. (22), 
which also indicates the good convergence of the ratio of exponential moments already at the 
charm mass scale in contrast with the = moments studied in Ref. |2^. We also expect 
that in the double ratio of moments used here, the radiative corrections to the gluon condensate 
effects (their expression for the two-point correlator is however available in the literature ||30[| ) 
are much smaller than in the individual moments , such that they will give a negligible effect in 
the estimate of the gluon condensate. This value obtained at the same level of a^-accuracy as 

30| , § from the ratio of 
claiming that the SVZ 
). Our value is also in 
and 



12 


14, 




P 





value |l[] has been underestimated by about a factor 2 (see c 

agreement with the more recent estimate (0.07 ± 0.01) GeV^ from the r-like sum rules 

I 



FESR m 



m e ' e 



1 hadrons. A more complete comparison of different determinations is 



done in Table 2. 



^For a recent review on the heavy quark masses, see e.g. 24, p5| 



Table 2: Predictions for the gluon condensate, for the different mass-sphttings (in units of MeV) 
and for the leptonic widths (in units of keV). We use as{Mz) = 0.118 ± 0.006 from LEP and 
r-decay. 



Observables 


Input 


Predictions 


Data / comments 


(a,G2)[GeV]^ x 10^ 


M^ - Mr,^ = 108 
^c.o.m _Mr = 440 

Xb 

Average 
Charmonium masses 

e~^e~ -^1 = 1 hadrons 

r-decay (axial) 

r-decay 
ALEPH 
CLEO 


10 ±4 
6.5 ±2.5 
7.5 ±2.5 

~ 4 

5.3 ± 1.2 

10 ±2 
9.2 ±3.4 

4± 1 

1 Q+5 

7± 1 
6.9 ±2.6 


This work 

Mass splittings 

SVZ-value 
g2 = 0-mom. 
g^-mom. 
exp. mom. [11, [T^ 
mom. 1 3 1 

ratio of mom. 33 
rJl/oK |34| 
r-like decay 

m 

data 
7 5+31 M 

2.0 ±3.8 [| 


a,(1.3 GeV) 




4p;+o-i8 

U-'±O-0.29 


a^iMz) ~ 0.124 ±0.012 


Xc(Pi) Xc(Po) 
^XrXPi) - ^Xc(Pf ) 


as from LEP/r-decay 
(ttsG^) average 


io.i;^i 

89l^,f 


11.1 (c.o.m) 

95 
50 


Mt - Mr,, 
^^■^Xb(P§) ^^^^ 

M^UP?) - Mr 
^x.(P|) - Mr 

^^^Xb(Pl) ^^^^ 


(asC^) average 
c.o.m. 


13+10 
63-29 

OJ+Sl 

475+Z^ 

^' '-'-64 

485ti 
500 ± 71 
492+^n 


order 
coef. as'. Coul. pot. 
400 

432 
453 
440 


Mr - 2mt 
Mr - Mr,, 

M^, - Mr 


{asG"^) average 


-906 
1.8 
93 

1800 


two-loop pole mass 

order 
coeff. as- Coul. pot. 


Tr^e+e- 




1.2 
0.16 


1.32 



Charmonium P-wave splittings 



The analysis of the different ratios of moments for the P-wave charmonium shows |]TT|-||I^ that 
they are optimized for: 



a - (0.6 ±0.1) GeV" 



aAa) 



0.41 



+0.11 
0.07 



l/x = 4.5 ± 1.5. (23) 



In the case of the Scalar Pq - axial Pf mass splitting, the different exact QCD coefficient 
corrections of the corresponding double ratio of moments read: 

= 1, ^ -(0.045;°:°^^), /SP^ ^ -(7.75l|f,). (24) 

Using the correlated values of the different parameters, one obtains the mass-splitting AMj^Q = 
Mp3 — Mpz ~ MeV, where we have used the experimental value Mps = 3.51 GeV. Adding 

the (gG^) dimension-six condensate effect,which is about -1.6% of the leading term in T^oi? one 
can finally deduce the prediction in Table 2, which is in excellent agreement with the data . One 
should remark that the previous predictions indicate that, for the method to reproduce correctly 
the mass-splittings of the 5* and P-wave charmonium states, one needs both larger values of 
and {asG"^) than the ones favoured in the early days of the sum rules. 

In the case of the Tensor P2^-axial Pf mass splitting, the different exact QCD corrections 
for the double ratio of the tensor over the axial meson moments read: 

A^^ = {0.989t'oZ) = (0.029;°:°?^) A^^ = (22.ll|^), (25) 

from which, one can deduce the prediction in Table 2, which is slightly higher than the data of 50 
MeV. This small discrepancy may be attributed to the unaccounted effects of the dimension-six 
condensate or/and to the (usual) increasing role of the continuum for state with higher spins. 
However, the prediction is quite satisfactory within our approximation. 



as from the - Pf axial mass splitting 



The corresponding double ratio of moments has the nice feature to be independent of the gluon 
condensate to leading order in and reads: 



^1 



A^3 (exact) = 0.014;°; 



008 



(26) 



The recent experimental value of the state denoted by hc{lP) in the PDG compilation p5 
coincides with the one of the center of mass energy: 



Mpi ~ Mi 



5Mp3 + 3Mp3 + Mpi 



11.1 MeV 



(27) 



as expected from the short range nature of the spin-spin force. We use the experimental value 
of the hc{lP) mass of 3526.1 MeV, and a naive exponential resummation of the higher order a.^ 
terms. Then, we deduce: 



1.3 GeV) ~ 0.64+°:?^ ±0.02 



as{Mz 



z) 0.127 ± 0.009 ± 0.002 ± 0.032 ± 0.030, 

(28) 

where The error is much bigger than the one from LEP and r decay data, but its value is perfectly 
consistent with the later. The theoretical error is mainly due to the uncertainty in Aq, while a naive 



resummation of the higher order terms leads the second error. However, though inaccurate, 
this value of as is interesting for an alternative derivation of this fundamental quantity at lower 
energies, which can serve for testing its g^-evolution until Mz- 

Reciprocally, using the value of ag from LEP and r-decay data as input, one can deduce the 
prediction of the center of mass (c.o.m) of the Pj states given in Table 2. 



T — r/5 mass splitting 



For the bottomium, the analysis of the ratios of moments for the S and P waves shows that they 
are optimized at the same value of a, namely [p!2[ : 

(29) 



a = (0.35 ± 0.05) GeV" 



which implies for 5 flavours: as{cr) — 0.32 ± 0.06. Using the conservative values of the two-loop 
6-quark pole mass: m?, ^ 4.2 — 4.7 GeV, one can deduce: 



1/x ^ 28 ± 7 , 

where one might (a priori) expects a good convergence of the 1/m expansion. 



(30) 



The splitting between the vector T{Sf) and the pseudoscalar rib{SQ) can be done in a similar way 
than the charmonium one. The double ratio of moments reads numerically: 



n 



VP 



Ml 
Ml 



(0.9995+°°o°2^ 



-o.ooosJ 



1 + a, 



(2.4 



0.7 

+1.4 



X 



10^3 + (0.03 ± 0.01)GeV-^(a,G2^ 



(31) 



where we have used the exact expressions of the QCD corrections. It leads to the mass splitting 
in Table 2. To this order of perturbation theory, this result is in the range of the potential model 
estimates [^G], |2^ , with the exception of the one in and ||2^ , where in the latter it has been 
shown that the square of the quark velocity f ^ correction can cause a large value of about 100 
MeV for the splitting. One should also notice that, to this approximation, the gluon condensate 
gives still the dominant effect at the fe-mass scale (0.2% of the leading order) compared to the one 
.08% from the a^-term. However, the 1/m series of the QCD as correction is badly convergent, 
showing that the static limit approximation is quantitavely inaccurate in the 6-channel. Therefore, 
one expects that the corresponding prediction of (13;^{q) MeV is a very crude estimate. In order 
to control the effect of the unknown higher order terms, it is legitimate to introduce into the sum 
rule, the coefficient of the Coulombic effect from the QCD potential as given by the a^-term in 
Eq. (12)[]. Therefore, we deduce the "improved" final estimate in Table 2: 



Mr - M,, ^ (63;^9) MeV, 
implying the possible observation of the eta^ from the Upsilon radiative decay. 



(32) 



T — Xb mass splittings and new estimate of the gluon condensate 

As the S and P wave ratios of moments are optimized at the same value of a, we can compare 
directly, with a good accuracy, the different P states with the T (Sf) one. As the coefficient of 
the corrections, after inserting the expression of cXmm, are comparable with the one from the 
Coulombic potential, we expect that the prediction of this splitting is more accurate than in the 

^In this case, the gluon condensate contribution is smaller than the Coulombic one as has been observed in jss) . 



case of the hyperfine. The different double ratios of moments read numerically for the values in 
Eqs (28)-(29): 
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(33) 



where V, S, A, T refer respectively to the T and to the different Xb states Pg , Pf , Using 
the value of the gluon condensate obtained previously, these sum rules lead to the mass-splittings 
in Table 2, which is in good agreement with the corresponding data, but definitely higher than 
the previous predictions of ||3^, where, among other effects, the values of as and of the gluon con- 
densate used there are too low. Reciprocally, one can use the data for re-extracting independently 
the value of the gluon condensate. As usually observed in the literature, the prediction is more 
accurate for the c.o.m., than for the individual mass. The corresponding numerical sum rule is: 



^c.o.m _ 



M- 



T 



(l.53l°i^) X 10-^ + (l.20l°i) X 10-2 + (0.28l°:°«)GeV-^(a.G2) 



which leads to: 



,G^) ^ (6.9 ± 2.5) X 10-2 GeV^ 



(34) 



(35) 



We expect that this result is more reliable than the one obtained from the — M^^ as the latter 
can be more affected by the non calculated next-next-to-leading perturbative radiative corrections 
than the former. An average of the two results from the ip — rjc and T — Xb mass splittings leads 
to: 



{asG^) ~ (7.5 ± 2.5) x 10"^ GeV^ 



(36) 



where we have retained the most precise error. This result can be compared with different fits of 
the heavy and light quark channels given in Table 2, and which range from 4 (SVZ) to 14 in units 
of 10~2 GeV^. The most recent estimate from e~^e~ ^1=1 hadrons data using r-like decay is 
(7 ± 1) X 10-2 GeV^ where one should also notice that the different post-SVZ estimates favour 
higher values of the gluon condensate. More accurate measurements of this quantity than the 
already available results from r-decay data |^ are needed for testing the previous phenomenological 
estimates from the sum rules. 



Toponium: illustration of the infinite mass limit 

Because in the alone case of the toponium, the 1/m-expansion is ideal, we have extended the 
previous analysis in this channel, though, we are aware that this application can be purely academic 
because of the eventual inexistence of the corresponding bound states. We use the top mass: 
rrit — (173 ± 14) GeV, obtained from the average of the CDF candidates (174 ± 16) GeV and of 



the electroweak data (169 ± 26) GeV as compiled by PDG [^. We shall work with the ratio of 
moments in the vector channel for determining the mass of the Sf state, and use with a good 
confidence the leading terms of the expressions given in Table 1. Using the sum rules in Eqs. (8) 
and (11) and the value of the minimum cj-^/^ ~ 20 GeV from Eq. (9), we deduce the result for the 
meson-quark mass gap given in Table 2. For the splittings, we use the sum rules in Eqs. (12) and 
(13), while, for the leptonic width, we use the sum rule in Eq. (14). Our results are summarized 
in Table 2. 



Conclusions 

We have used new double ratios of exponential sum rules for directly extracting the mass-splittings 
of different heavy quarkonia states. Therefore, we have obtained from — M^^ and M^^^ — Mj- 
a more precise estimate of the value of the gluon condensate given in Eq. (34). We have also 
used M^^(pi) — M^^(p3) for an alternative extraction of at low energy (see Eq. (28)), with a 
value consistent with the one from LEP and r-decay. Our numerical results are summarized in 
Table 2, where a comparison with different estimates and experimental data is done. We have 
not extended this analysis to the D-waves as these states are described with operators of higher 
dimensions, where, in general, the role of the continuum is relatively important and can obtruct 
the extraction of the mass-splittings. We have also attempted to connect the sum rules and the 
potential model approaches, using a 1/m-expansion. We found, that the Coulombic corrections, 
which are quite well understood in QCD, agree, in general, in the two approaches, expect in the 
radiative corrections of the hyperfine splitting which requires the knowledge of the next-next-to- 
leading a^-corrections. Relativistic corrections due to finite value of the quark mass have been 
included in our analysis. However, the coefficients of the gluon condensate disagree in the two 
approaches, which may be related to the difficulty encountered by Bell-Bertlmann in finding a 
bridge between a field theory a la SVZ and potential models. 
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